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ABSTRACT 

The intrinsic shape of galaxy clusters can be obtained through a combination of X-ray and Sunyaev- 
Zeldovich effect observations once cosmological parameters are assumed to be known. In this paper we discuss 
the feasibility of modelling galaxy clusters as either prolate or oblate ellipsoids. We analyze the intra-cluster 
medium distribution for a sample of 25 X-ray selected clusters, with measured Sunyaev-Zeldovich temperature 
decrements. A mixed population of prolate and oblate ellipsoids of revolution fits the data well, with prolate 
shapes preferred on a ~ 60 — 76% basis. We observe an excess of clusters nearly aligned along the line of 
sight, with respect to what is expected from a randomly oriented cluster population, which might imply the 
presence of a selection bias in our sample. We also find signs that a more general triaxial morphology might 
better describe the morphology of galaxy clusters. Additional constraints from gravitational lensing could dis- 
entangle the degeneracy between an ellipsoidal and a triaxial morphology, and could also allow an unbiased 
determination of the Hubble constant. 

Subject headings: Galaxies: clusters: general -X-rays: galaxies: clusters - cosmology: observations - distance 
scale - gravitational lensing - cosmic microwave background 
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1. INTRODUCTION 

Measuring the three-dimensional shape of any class of as- 
tronomical object is today still an almost unsolved task which 
has challenged astronomers for decades. The knowledge of 
the shape of clusters of galaxies in particular has fundamen- 
tal cosmological implications. In a cold dark matter scenario, 
smaller clumps aggregate hierarchically along large-scale per- 
turbations, which are typically highly anisotropic at their first 
collapse. Accretion of material might hence preferentially oc- 
cur along the filamentary structure within which larger matte r 
aggregates (i.e. galaxy clusters) are embedded feWesll l1994l) . 
This could explain the observed tendency of clusters to be 
aligned with their nearest neighbor (Carter & Metcalfe 1980). 
Further processes, such as virialization of the matter aggre- 
gates or self-interaction mechanisms for dark matter, might 
on the other hand tend to make such systems more spheri- 
cal. Dissipation and gas cooling could result in an average in- 
crease of the axis ratios of dark matter halos by 20 — 40% 
in the inner regions, with effects persisting almost to the virial 
radius (Kazantzidis et al. 2004). The intrinsic shape of struc- 
tures therefore contains evidence of the formation history of 
the large-scale structure and abo ut the nature and mechan isms 
of interaction of dark matter tfhonis et £0119911 120041 and 
references therein). 

The first attempts to determine the three dimen- 
sional morphology of galaxies with only photometric data 
(jBinnev & Merrifield 1998, and references therein) led to the 
reaUzation that since objects might be either triaxial or build 
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up a mixture of oblate and prolate shapes, a purely geometri- 
cal and stat istical analysis will not produce a unique solution 
('Bing2eli'1980). 

First results came from the orig inal work oflHubbld {^3), 
who first determined the relative frequencies with which 
galaxies of a given intrinsic ellipticity, oriented at random, 
will be observed as having various apparent projected elliptic- 
ities. Several studies have subsequently generalized Hubble's 
work ( Binnev & Merrifield 1998, and references therein), 
making the statistical approach of inversion of the distribu- 
tion of the apparent shapes of galaxies a widely used method. 
This inversion method has several drawbacks, though. First 
of all, it is not unambiguous since it requires strong assump- 
tions on the distribution of the intrinsic shapes of the struc- 
tures. The distribution of the apparent flattening is moreover 
rather insensitive to modifications in the model assumed for 
the intrinsic shapes iBinggeIiir9 80): different distributions of 
oblate or prolate spheroids, or generally triaxial ellipsoids, 
can consequently be consistent with the observed distribu- 
tion of axial ratios. Simplifying hypotheses, such as requir- 
ing structur es to be oblate, can he nce be strongly mislead- 
ing ( Bertola^^aEiccio^l^^llllingworth 1977); for a more 
correct approach, prolate and triaxial models should also be 
considered. 

While the first investigation by Hubble ( 1926) was done on 
galaxies, several following studies have faced the inversion 
method for different classes of ast ronomical objects such as 
poor groups or clusters of galaxies ( Binnev & de Vaucouleurs] 
1981; Thakur & Chaki-abortv 2001; Fasano & Vio 19911 
Noerdlinger 1979; Alam & Rvden 2002; BinggeM 198(J 
■Rvden 1992, .1996; Phonis et al. 2004) . With the exception 
of disc galaxies, prolate-hke sh apes appear to dorn inate all 
cosmic structure on a large scale dPlionis et a lJ2004ft . 

In this paper, we are interested in obtaining the three- 
dimensional shape of clusters of galaxies. So far, studies on 
clusters have been limited to the statistic al invers ion of the ap- 
parent distribution of axial ratios. Plion is et al J il99 1 ) found 
that the distribution of apparent axial ratios of 397 Abell clus- 
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ters, selected from the Lick map of galaxies, is inconsis- 
tent with a pop ulation of oblate sphero ids. Similar results 
are reported by iBasilakos et al.l J2000h and Ide Theiie et al.l 
([1995), who find that a population composed by purely pro- 
late spheroids is more consistent with, respectively, a sample 
of 903 clusters from the APM catalog and with a set of 99 
low redshift Abell clusters, better than a purely oblate pop- 
ulation. While the above studies were performed uniquely 
on optical data, Coorav (2000) used the X-ray i sopho tal axial 
ratios of a sample of 58 clusters by Moh r et aP J 1 9951) : a pop- 
ulation of purely oblate ellipsoids is again rule d out. A quite 
different approach was followed in Paz et al] where 
observational data from 2PIGG and SDSS group catalogues 
were compared to projected shapes of numerically simulated 
haloes. Haloes turned out to be preferentially prolate, with 
more massive groups tending to show more elongated shapes. 

The above listed results are based on the controversial as- 
sumption that objects are randomly oriented in space. Fur- 
thermore, tests for possible intrinsic triaxial or mixed popu- 
lations have to date been far from conclusive. Several works 
in the past few years have discussed how multi-wavelength 
observations of clusters of galaxies can be used to uncover 
the 3- D structure of clusters of gala xies ijZaroubi et alJll99^ 
1200 U Reblinskv 2000; Fox & PenI QooIk I n a com pan- 
ion paper (De Filippis et al. 2005) (hereafter. iPaper ll) . we 
have considered the theoretical capability of combined X- 
ray, Sunyaev-Zeldovich effect (SZE) and gravitational lens- 
ing observations to determine the 3-D morphology of relaxed 
galaxy clusters and, jointly, to constrain cosmological param- 
eters. Due to a lack of gravitational lensing data, we have then 
used X-ray and SZE measurements of a sample of 25 clusters 
to constrain their triaxial structure by assuming the concor- 
dance cosmological ACDM model, which is determined to an 
unprecedented accuracy thanks to an impressive body of ev- 
idence coming from several cosmological tests iWang etal] 
EOOO; Tegmarket al. 2004). 

In this paper, we test the accuracy to which ellipsoids of 
revolution can des cribe the sample of clusters of galaxies we 
studied in iPaper S Following a similar empirical approach, 
we use combined X-ray and SZE observations to directly 
probe the 3-D structure of clusters in the sample, without 
any assumptions on their statistical properties. The paper is 
organized as follows. In § |2l we discuss how the shape of 
an ellipsoid of revolution can be obtained from combined X- 
ray and SZE measurements. In §|3lwe introduce our clus- 
ter sample and derive the intrinsic distributions of the axial 
ratios and inclination angles for the cluster sample, together 
wi th a com parison between these results and those obtained 
in IPaper II assuming a triaxial morphology aligned along the 
line of sight. In §|3we perform a statistical inversion of the 
observed distribution of projected axial ratios. In § |6l some 
systematics are considered. Further information that can be 
obtained with additional constraints from gravitational lens- 
ing observations is discussed in § [S] § is devoted to con- 
clusions and final considerations. In AppendixlAl we provide 
some mathematical details on the two-dimensional projection 
of an ellipsoid. 

2. ELLIPSOIDS OF REVOLUTION 

In this paper we assume that galaxy clusters have an ax- 
isymmetric gas density distribution. The intrinsic distribution 
of the intra-cluster medium (ICM), assumed to be constant 
on similar, concentric spheroidal ellipsoids, is characterized 
by the ratio of the major to minor axes eint(> 1), and by 



the inclination angle i, between the line of sight (LOS) of the 
observer and the polar axis. In an intrinsic orthogonal coor- 
dinate system centred on the cluster's barycentre and whose 
coordinates are aligned with its principal axes, a spheroidal 
ICM profile can be described by only one radial variable C,, 
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If we take the axis of symmetry aligned with the third co- 
ordinate, ei = {eint,eint, 1} for a prolate model and = 
{1,1, eint} for an oblate model. 

For an intrinsic 3-D spheroidal morphology, the two- 
dimensional observed quantities, projected on the plane of 
the sky (POS), are constant on concentric ellipses. The ax- 
ial ratio of the major to minor axes of the observed projected 
isophotes, eproj, can be expressed a s a function of the cluster 
intrinsic parameters dFabricant et aU1984) . 
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In lPaper fl we discussed how the elongation along the LOS 
of a cluster of galaxies, elos, can be determined by combin- 
ing X-ray and SZE observations if a reliable estimate of the 
Hubble constant is provided. For a spheroid inclined with re- 
spect to the LOS, as we assume clusters to be shaped in this 
paper, glos represents the ratio of the size of the spheroid 
along the LOS, measured through its barycentre, to the size of 
the major axis projected in the POS. A value of glos greater 
than one corresponds to a cluster elongated along the LOS, 
and vice versa for glos < 1- For the cluster plasma density 
we assume an isothermal /3-distribution: 

/ a2 \ -3/3/2 

ne(C) CX 1 + \ , (3) 



where is the core radius and (3 the slope parameter. For 
such a distribution, glos can then be written as (Paper L) : 
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r[3/3/2- 1/2] 
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where z is the cluster redshift, D^. the angular diameter dis- 
tance to the cluster, ATg the central SZE temperature decre- 
the projected angular major core radius, 



ment. 



^c.proj 



the temperature of t he ICM, kn the Boltzmann constant, 
TcMB = 2.728°K dFixsen et all Fl 996.) the temperature of 
the cosmic microwave background radiation, itt the Thomp- 
son cross section, mo the electron mass, /imp the mean par- 
ticle mass of the gas, c the speed of light in vacuum, Ae 
the X-ray cooling function of the ICM in the cluster rest 
frame, Sxo the central X-ray surface brightness and /(i/. To) 
accounts for frequency shift and relativistic corrections. In 
a flat Friedmann-Lemaitre-Robertson- Walker model of uni- 
verse filled with dust, ^Imo, and with a non null cosmological 
constant, the angular diameter distance to a source at a red- 
shift z is (.Sereno et al. 200'll) 

- ^ W\ (5) 
Hol + z Jo ^S1mo(1 + + 1 - ^mo 
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We assume a flat ACDM model with a Hubble constant Hq ~ 
70tt km s-^Mpc-i and nuo = 0.30±0.04 ( Tegmai-k et al. 
120041) . As well known, if the Hubble constant is not assumed 
from independent observational constraints, then the length of 
a cluster along the LOS can be known from SZE and X-ray 
observations only up to a multiplicative factor (Coorav 199 8). 
When we use experimental data pro vided with asymm etric 
uncertainties, connections as given bv lD'Agostinil J2004 1) are 
applied to obtain unbiased estimates of mean and standard 
deviation. 

In terms of en mi w e can also write (iDe Filippis et al.l2005t 



iDonahiie etal.i2003l) 



glos 



proj 
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prolate case 



oblate case 
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For clusters aligned along the LOS (i.e. i — 0) both prolate 
and oblate distributions exhibit circular isophotes and glos 
can be expressed as a function of the cluster intrinsic elliptic- 
ity as follows: 



eint prolate case 
oblate case. 

eint 



(7) 



If we assume the cluster to be either prolate or oblate 
Eqs. enable us to determine the intrinsic geometry of 

the cluster. In the prolate case, we obtain 

^eLoseproj; 



eint = 



(8) 



: ± arccos 



while in the oblate case: 

en 



eint 



CLOS 



: ± arccos 
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We are dealing with observed projected quantities which are 
invariant under a reflection through the POS. Hence we are not 
able to determine which extremity of the cluster is pointing 
towards the observer. 

3. INTRINSIC DISTRIBUTIONS THROUGH COMBINED X-RAY AND 
SZE MEASUREMENTS 

In this section, we apply the formalism described in §|2]to 
compute the three-dimensional intrinsic elongation and incli- 
nation for the same cluster sample used in Paper I, for which 
a combined X-ray and SZE data set was available. Th e sam- 
ple co nsists of 25 clusters of g a laxies , 18 from Rees e et al.l 
(I2002h and 7 from 'Mason et al? f2001'). Details on the sam- 
ple an d on the selection criteria can be found in Mason et al. 
boOlh : iReese et alJ (2002); De Filippis et al. (2005). Basic 
data for our 25 clusters, including previously published red- 
shift, plasma temperatur e and SZE decrement information are 
presented in Table 1 of ^ Paper li In iPaper I. we also mea- 
sured projected axial ratios, Cproj, for all clusters in the sample 
through an X-ray morphological analysis based on high res- 
olution Chandra and XMM- Newton observations; measured 
values are listed in Table 2 of iPaper ll 



Through a combined analysis of both X-ray and SZE ob- 
servations, we can infer new information about intrinsic elon- 
gation and inclination of the clusters in the sample. We 
then m odel these distributions with a kernel empirical esti- 
mator ( Vio et al. 1994; Ryden 1996). Given a sample of N 
measured quantities r/i, the kernel estimator of the distribu- 
tion is a continuous function given by: 



1 ^ 



(12) 



where K is the kernel function. When densities take values in 
the range < 77 < rymax, the kernel can be written as: 



Gau 



h J V h 



where ifcau is a Gaussian kernel given by: 



1 



(14) 



When using the kernel in Eq. il3\ . Gaussian tails extending 
to negative values or to values greater than r^max, are folded 
back between and ?7niax- The kernel width h can be approx- 
imated by /i ~ 0.9AN~'^-^, where A is the smallest quan- 
tity between the standard deviati on of the sample and the in- 
terquartile range divided by 1.34 d Vio et alJl994ft . 

In order to account for the effect of the finite sample size, 
we derive the confidence levels from a bootstrap re-sampling 
of the data set. From the distribution obtained with the orig- 
inal N data points, we draw with replacement a new set of 
N data points which are then used to create a new bootstrap 
estimate of the kernel estimator, iV^. Confidence intervals can 
then be assigned to each value of 77, by finding the values of 
that lie above (below) the required upper (lower) confi- 
dence limit. Errors in the measured value of 7/ induce a fur- 
ther uncertainty in the estimated distribution. This effect can 
be modelled with a kernel width given by: 



(15) 



which changes for each data point dRvdenll996h . 

In what follows, we use the inverse of the intrinsic and pro- 
jected axial ratios: qi^t = 1/eint and qpioj = 1/eproj- Such a 
formalism allows axial ratios to be limited to the finite range 
(0, 1]. The histogram in Fig. ^ represents the observed dis- 
tribution of the projected axial ratios. Confidence limits are 
built up with 10'^ bootstrap re-samplings. Estimates of the 
observed projected ellipticity from both Chandra and XMM- 
Newton satellites are so accurate that for the distribution of 
observed axial ratios ~ h. The measured mean ellipticity 
is 9pioj = 0.80±0.02; the value of the median is instead 0.81. 
Reported uncertainties represent the standard deviation of the 
mean. These values are in good agreement with the average 
projected axial ratio expected if dark matter halos of clusters 
do have mass distributions as revealed by numerical simula- 
tions (Wang & Fan 2004). The solid and the dashed lines in 
Fig- n show the non-parametric kernel estimate and the 2-0- 
confidence range, respectively, using the bootstrap procedure. 

3.1. Axial ratios 
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TABLE 1 

Intrinsic Parameters of Sample Clusters 









<3int 


i 


9int 




Cluster 


z 


Shape 


(ellips.) 


(deg) 


(triax.) 




MS 1137.5+6625 


0.784 


pro 


0.49 ± 0.21 


16 ±9 


0.56 ± 


21 


MS 0451.6-0305 


0.550 


obi 


0.71 ±0.16 


58 ±43 


0.76 ± 


06 


- 


- 


pro 


0.63 ±0.14 


56 ± 13 


- 




CI 0016+1609" 


0.546 


pro 


0.55 ±0.13 


26 ± 10 


0.67 ± 


14 


RXJ1347.5-1145 " 


0.451 


obi 


0.68 ±0.16 


84 ± 117 


0.69 ± 


09 


- 


- 


pro 


0.48 ±0.11 


35 ± 12 


- 




A 370 


0.374 


pro 


0.35 ±0.13 


27 ± 12 


0.55 ± 


15 


MS 1358.4+6245 " 


0.327 


obi 


0.55 ± 0.14 


52 ± 8 


0.71 ± 


06 


A 1995 


0.322 


obi 


0.75 ±0.18 


64± 35 


0.81 ± 


09 


- 


- 


pro 


0.70 ±0.16 


46 ±27 


- 




A611 


0.288 


pro 


0.73 ± 0.26 


35 ± 32 


0.84 ± 


19 


A 697 


0.282 


pro 


0.44 ±0.12 


26 ± 10 


0.60 ± 


14 


A 1835 •■■ 


0.252 


pro 


0.56 ± 0.08 


29 ±7 


0.70 ± 


09 


A 2261 " 


0.224 


pro 


0.63 ±0.16 


10 ±6 


0.65 ± 


16 


A 773 


0.216 


pro 


0.32 ± 0.09 


14 ±5 


0.40 ± 


11 


A 2163" 


0.202 


pro 


0.62 ±0.13 


28 ± 10 


0.73 ± 


12 


A 520 


0.202 


obi 


0.54 ± 0.20 


23 ± 10 


0.56 ± 


08 


A 1689" 


0.183 


obi 


0.85 ±0.12 


65 ± 44 


0.88 ± 


06 


- 


- 


pro 


0.79 ±0.11 


46 ± 23 


- 




A 665 


0.182 


obi 


0.58 ± 0.25 


47 ± 13 


0.71 ± 


15 


A 2218 


0.171 


pro 


0.54 ±0.18 


22± 11 


0.64 ± 


18 


A 1413" 


0.142 


obi 


0.64 ±0.17 


73 ± 33 


0.68 ± 


09 






pro 


0.49 ±0.13 


37 ± 15 






A 2142" 


0.091 


obi 


0.63 ±0.13 


79 ±42 


0.65 ± 


07 






pro 


0.43 ± 0.09 


34 ± 10 






A 478 " 


0.088 


pro 


0.34 ± 0.14 


23 ± 12 


0.51 ± 


17 


A 1651 " 


0.084 


pro 


0.31 ±0.16 


12 ±7 


0.37 ±0 


19 


A 401 


0.074 


pro 


0.46 ± 0.08 


25 ± 6 


0.60 ± 


10 


A 399 


0.072 


obi 


0.50 ±0.14 


38 ± 5 


0.58 ± 


04 


A 2256 


0.058 


pro 


0.54 ±0.14 


34 ± 14 


0.72 ± 


11 


A 1656 


0.023 


pro 


0.70 ± 0.28 


33 ± 30 


0.81 ± 


21 



Note. — For ambiguous clusters, both prolate and oblate solutions are listed, z 
is the redshift. qint (ellips.) the cluster maximum axis ratio under an ellipsoidal as- 
sumption, i the inclination angle, gint (triax.) the cluster maximum axis ratio under a 
triaxial assumption. 

^Cooling flow clusters. 
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Fig . 1 . — Normalized distribution of apparent axial ratios. Solid and daslied 
lines sliow the non-parametric kernel estimate and the 2-a confidence range 
using a bootstrap procedure. 



In iPaper J . the approach used to calculate the elongation 
along the LOS is discussed in detail; estimated values of eLos 
are listed there in Table 4. Once eLos and epioj are known, 
both systems Eqs. (|8j |9} and Eqs. dlOl II W for prolate and 
oblate ellipsoids, respectively, can be solved. Given each set 
of observational constraints, a cluster can admit either only 
one solution (prolate or oblate) or both. By solving the equa- 
tion systems, we find that of the 25 clusters in our sample, 15 




0.2 0.4 0.6 0.8 1 

qint 



Fig. 2. — Distribution of intrinsic axial ratios. Histogram and thick lines 
are computed assuming a prolate solution for all ambiguous clusters. Thick- 
solid and thick-dashed lines show the non-parametric kernel estimate and the 
2-cr confidence range using a bootstrap procedure. The thin line is the kernel 
parametric estimate assuming the oblate solution for all ambiguous clusters. 
Histogram and all distributions are normalized to unity. 



are compatible only with a prolate morphology, 4 only with 
an oblate one, while the remaining 6 clusters (hereafter the 
"ambiguous clusters") are compatible with both the prolate 
and the oblate assumption. The inferred intrinsic parameters 
are listed in Table ^ For comparison, in the last column of 
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Fig. 3. — Distribution of inclination angles. Histogram and lines are the 
same as in Fig.l2l 




Fig. 4. — Orientation angle versus the intrinsic ellipticity for all clusters in 
our sample. Boxes and triangles represent prolate and oblate clusters, respec- 
tively. Both solutions for ambiguous clusters are plotted in gray. 



Table [2 we list the values of q-mt as computed in [Paper 3 as- 
suming clusters as triaxial structures with null inclination. 

To analyze the inferred distribution and to model it with 
the kernel estimator, we first assume that all of the ambigu- 
ous clusters are prolate (Case I) and then that they are all 
oblate (Case II). These are both extreme cases, but the in- 
trinsic distribution turns out not to be very sensitive to this 
assumption. In Fig.|2]we plot the distribution of the intrin- 
sic axial ratios. Also in this case, confidence ranges have 
been determined using the bootstrap procedure with 10^ re- 
samplings. It can be seen that the kernel estimate does not 
change significantly if all ambiguous clusters are assumed to 
be all prolate or oblate (solid thick and thin line, respectively). 
The distribution obtained in Case II (thin line) is well within 
the 2-CT confidence level of the one obtained in Case I (thick- 
dashed lines). In Case I (II), the mean intrinsic axial ratio is 
0.54±0.03 (0.56±0.03), while the median is 0.54 (0.56). 

3.2. Inclination angles 

The polar axis of clusters is usually assumed to be randomly 
oriented with respect to the line of sight. Under this assump- 
tion a fraction sin i di of all clusters would have their sym- 
metry axes directed within di of angle i with respect to the 
line of sight (Binnev & Merrifield 1998). The mean value of 
such a distribution would be (i) random ^ 57 . We instead ob- 
serve a slightly different picture. The measured distribution 
of the inclination angles is shown in Fig.|3] The mean incli- 
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Fig. 5. — 3-(j confidence bands of the intrinsic axis ratio distribution, 
as determined through combined X-ray and SZE measurements with the as- 
sumption of an inclined ellipsoidal morphology (dark shaded region. Case I), 
and as determined under the assumption of triaxial morphology aligned along 
the LOS (light shaded region). Confidence bands are found by bootstrap re- 
sampling. 



nation angle is 32±3 deg (37±4deg) for Case I (II), whereas 
the median is 29 deg for both cases. Observed clusters show 
an excess of clusters with their polar axes closely aligned with 
the line of sight, and a strong deficit of clusters with high val- 
ues of i. The distribution plotted in Fig.lsjhas been obtained 
without considering the error Ai when calculating the kernel 
widths. Using instead of hi causes a further broadening of 
the distribution; both the excess of small values of i and the 
discrepancy with the assumption of a random distribution still 
hold. 

Fig- S shows a plot of the orientation angle versus the in- 
trinsic ellipticity for all clusters in our sample. The most 
elongated clusters (small values of qi^t) turn out to be pro- 
late nearly aligned with the line of sight. This is most prob- 
ably due to a bias in the selection of the clusters in the sam- 
ple. Clusters with low X-ray luminosity, which are stretched 
along the line of sight are preferentially included in a surface 
brightness limited survey. Both the ellipticity and the inclina- 
tion inferred for ambiguous clusters are quite sensitive to the 
assumed shape. The oblate solution implies a more spherical 
geometry (larger values of gint) and a more inclined polar axis 
than does the prolate solution. In fact, ambiguous clusters can 
have inclination angles > 60 deg only if they are oblate. 

3.3. Rotational versus triaxial spheroids 

While in lPaper 3 we assumed galaxy clusters to be triaxial 
ellipsoids ahgned along the LOS, in this paper we have re- 
laxed the assumption of the null inclination, allowing clusters 
to be inclined along the LOS, but we have put a slightly tighter 
constraint on their morphology, forcing them to be ellipsoids 
of revolution. 

We observe a tight concordance between the results ob- 
tained in the two analyses. Estimated values of gint computed 
under the two assumptions are in agreement in most cases (see 
Table^. The case of a cluster compatible only with the pro- 
late (oblate) shape corresponds to a triaxial cluster with the 
major (minor) axis oriented along the line of sight. Ambigu- 
ous clusters correspond to triaxial shapes having the middle 
axis directed along the line of sight. The two distributions of 
the intrinsic axis ratio as obtained either in the present analysis 
or in Paper I are represented in Fig. |5l (dark and light shaded 
regions, respectively). As a general trend, we observe that 
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Fig. 6. — Deprojected distribution of intrinsic axis ratios assuming tliat all 
clusters are randomly oriented prolate spheroids (solid line). Short-dashed 
and long-dashed lines show the 1 and 2-a confidence limits, respectively, 
estimated using a bootstrap re-sampling technique. 
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Fig. 8. — 3-(7 confidence bands of the intrinsic axial ratios, as determined 
from combined X-ray and SZE measurements (dark shaded region), and of 
the deprojected distribution of observed axial ratios, assuming that clusters 
are randomly oriented prolate spheroids (light shaded region). Confidence 
bands are found by bootstrap re-sampling. 




Fig. 7. — Deprojected distribution of intrinsic axial ratios assuming that 
all clusters are randomly oriented oblate spheroids (solid line). Short-dashed 
and long-dashed lines show the 1 and 2-u confidence limits, respectively, 
estimated using a bootstrap re-sampling technique. 



galaxy clusters turn out to be slightly more elongated when 
modelled as ellipsoids of revolution than in the case of triax- 
ial ellipsoids (showing an average ratio of the minor to the 
major axis of ^ 0.63±0.03). This effect is due to the capabil- 
ity of triaxial shapes to fit observed data without the need of 
extreme axis ratios. 

4. DEPROJECTED DISTRIBUTION 

In this section, we compute the distribution of the intrinsic 
axis ratio with a third method. The observed projected ax- 
ial ratio distribution, /((/proj)- can in fact be inverted under 
appropriate assumptions in order to obtain the intrinsic dis- 
tribution of the axis ratios. If clusters were all randomly ori- 
ented prolate ellipsoids, the intrinsic distribution f^-pm could 
be written as (Rvden 1996) 



^Pro(gint) 



2^1 



9L 



dq 



(16) 

If, on the contrary, all systems are assumed to be randomly 
oriented oblate ellipsoids, then the intrinsic distribution can 



be written as: 



A^Ob^gint) 




dq 



(17) 

To take into account the effect of the finite sample size, con- 
fidence levels for iVpio and iVobi have been derived using a 
bootstrap re-sampling procedure, see § |3l The bootstrap es- 
timates of / are then inverted to compute estimates of the 
intrinsic distributions Npm and iVobi- Confidence intervals 
are assigned to each value of gint by finding values of Np^o 
and Noh\ that lie above (below) some upper (lower) confi- 
dence limits. Accurate confidence limits were built using 
10'^ bootstrap re-samplings. The resulting intrinsic axial ra- 
tio distributions are plotted in Figs. |S] and these were ob- 
tained assuming purely prolate and oblate morphologies, re- 
spectively. Negative values in the inverted distri butions are 
evidence against the underlying hypotheses (Rvden 
ICoorav 2000). The assumption that all clusters are randomly 
oriented prolate ellipsoids turns out to be consistent with the 
observed distribution of projected axial ratios, see Fig.|6l On 
the other hand the hypothesis that all clusters are purely oblate 
is discarded at the l-a confidence limit since the resulting 
intrinsic distribution is negative for high values of qint, see 
Fig-El However, at the 2-cr confidence level, such an hypoth- 
esis cannot be rejected, given the observed ellipticities of our 
sample clusters. 

We now compare the intrinsic distributions of axis ratios 
obtained using our deprojection technique through X-ray and 
SZE observations, with the above inverted distribution. We 
have taken into account the 3- a confidence levels of the intrin- 
sic axis ratio distribution obtained with combined X-ray/SZE 
method, estimated independently for Case I and II in § [S] 
We have then obtained the concordance region by holding all 
points compatible with at least one of the cases. The resulting 
region is plotted in Fig.|8]in dark gray. The light gray region 
shows the inverted distribution of observed axial ratios, as- 
suming that clusters are randomly oriented prolate spheroids. 
With respect to our X-ray/SZE technique, the deprojection 
method leads to a deficit of nearly spherical structures (low 
values of qiat) and to an overabundance of highly elongated 
clusters. The same discrepancy is observed if the distribu- 
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tion of observed axial ratios is inverted assuming that clusters 
are randomly oriented oblate spheroids. Both assumptions 
that clusters are either all randomly oriented prolate or oblate 
spheroids do not provide suitable explanation for the intrinsic 
distribution of our sample. 

5. GRAVITATIONAL LENSING CONSTRAINTS 



As already discussed in lPaper fl. a third independent obser- 
vational input can both break the degeneracy between oblate 
and prolate models and determine cosmological parameters. 
Such an input can be provided by observations of strong lens- 
ing events. Whereas the ICM distribution traces the gravita- 
tional potential, gravitational lensing directly maps the cluster 
total mass. We remark that since we are dealing with an el- 
lipsoidal ICM distribution, the isodensity contours will not 
in general be ellipsoidal, this being approximately true only 
for a small value of the ellipticity t iFox & Penll2002^ . Obser- 
vations of strong lensing events can put accurate constraints 
on the convergence, i.e. the dimensionless projected surface 
density of the total mass distribution of the lens. Assuming 
hydrostatic equilibrium and a /j-model for the ICM , the cen- 
tral value of the convergence reads dPe Filippis et alJI2005L 
Eq. (C7)) : 



Cint 



1 



if prolate 



1 



c.proj 



^ if Oblate 

Bint 

(18) 

where D^s and Dg are the angular diameter distances from the 
deflecting cluster to the lensed background source and from 
the observer to the source, respectively. 

The equation for the convergence can be coupled to those 
discussed in §|3to determine if the cluster is oblate or prolate, 
and to estimate one cosmological parameter, in particular Ho, 
together with the shape parameters i and ejnt- Additional mul- 
tiple image syst ems allo w us to estimate further cosmological 
parameters (see lPaper^ . 

6. SYSTEMATICS 

Several systematic uncertainties that depend on the physical 
state of the ICM could affect our determination of the intrinsic 
structure of a galaxy cluster. Let us review some effects. 

6.1. Cooling flows and temperature gradients 

Departures from isothermality such as a non constant tem- 
perature profile or temperature sub-structures mainly affect 
our analysis through the SZE, since AT cx J n^Tdl. 

ICM temperature profiles have been accurately measured 
thanks to new g eneration X-ray satellites for several nearby 
relaxed clusters (iPiffaretti et alJ2005tl\^lin in et al.'' 2005bh . 
although, due to restricted field of view or background lev- 
els, in most cases such measurements are only possible out to 
one-third/half of the virial radius. Results are still controver- 
sial, but some recent investigations seemed to find a similar 
behaviour in the ICM temperature distribution, with a broad 
peak followed by a decrease at larger radii ( Piffaretti et al. 
l2005t 'Vikhlinin et al.' "2005aV Cooling processes also al- 
ter the temperatur e profiles with sharp drops within cluster 
cores J Allen et al.l 12001). Similar measurements for large 
samples of more distant clusters are instead still lacking. 

Since the SZE depends essentially on the pressure profile of 
the cluster, the presence of a cooling flow can affect our esti- 
mate of the cluster extent along the line of sight. Our assump- 
tion of isothermaUty in the X-ray analysis generally leads to 



an overestimate of the temperature and (to a lesser extent) the 
pressure in the cool region. This in turn leads to an overesti- 
mate of the expected SZE signal. The observed temperature 
decrement is thus misinterpreted, leading to an underestimate 
of the LOS extent. 

A refined view of cooling processes should account for 
changes in density profiles as well as the temperature profile, 
with the effect due to decrease of gas temperature moderated 
by the sudden increase of the gas density. Previous studies 
addressed the effect of departures from isothermality when 
using SZE along with X-ray observations. Given our differ- 
ent approach, we interpret as an over (under)-estimation of 
the LOS elongation what has been previously seen as an un- 
der (over)-estimation of the Hubble constant. Based on the- 
oretical models for radiative cooling, it was found that, even 
after excluding ~ 80% of the cooling-flow region, an error of 
~ 10% on the elongation could affect analyses restricted to 
radii of 1 . 5rc (iMaiumdar & Nath>.2000.) . 

In order to check how much the presence of a central cool- 
ing region can affect our analysis, we re-analyzed three mas- 
sive cooling flow clusters in our sample (RX J1347.5-1 145, 
A 1835 and A 478) using the most recent Chandra obser- 
vations available. We performed morphological and spectral 
analyses both using the whole cluster extent and excluding 
the central cool region. Details on the data reduction and on 
the morphological analysis are given in Paper I. Removing 
the central region leads to an increase of both values of 9c 
(« 60%) and j3 (> 6%), when fitting a /3-model to the clus- 
ter surface brightness elliptical profile. Estimates of central 
amplitude are strongly affected too. Average isothermal spec- 
tral temperatures are larger by ^ 20% respect to when the 
values are averaged over the whole cluster extent. The above 
effects leads to an underestimation of the cluster elongation 
along the line of sight of ^ 15%, in agreement with pre vious 
estimates fl nagaki et al.ll995H Maiumdar & Nath''2000). 

The effect of a central decrement is opposite to that of 
a negative gradient at large radii in the temperature profile 
when determinin g the c luster structure as suming isothermal- 
ity (gnagaki et all ll995t IVoshikawa et alJ 119981) . If the gas 
temperature in the outer region is lower than that in the in- 
ner region, an overestimation of the LOS extent as large as 
- 20% is possible (llnagaki et alJl995l) . 

Actually, systematic errors due to departures from isother- 
mality are overwhelmed by large observational uncertainties 
in both the actual temperature profile and the SZE measure- 
ment. Results turn out to be quite insensitive to details of the 
gas modelling and the isoth ermal /3-model ca n still provide a 
reasonably accurate model iBonamente et al.1.200 5). Effects 
due to either central drops or temperature gradients at large 
radii could partially compensate each other and could be com- 
bined in a systematic deviation as large as 15%, whose sign 
depends on which effect prevails. 

6.2. Density clumping 

Small-scale density fluctuations on X-ray measurements 
arising from accretion events and major mergers could intro- 
duce a further bias. Clumping in the gas density distribution 
causes an enhancement of the X-ray brightness by a factor 



with respect to a uniform smooth atmo- 



sphere, while SZE measurements are not affected. This leads 
to an underestimate of the elongation along the line of sight 

for Cn > 1 (CLOS oc C~i) . 

Currently, there is no observatio nal evidence of si gnif- 
icant clumping in galaxy clusters (iReese et alJ l2002h . al- 
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though nu merical hydr o-simulations show that Cn — 1.34 
ifMathiesen et al.lll999h . More detailed simulations and ob- 
servations are required to estimate the extent and the effect 
that dumpiness could induce on our analysis. 

6.3. Overall departure 

In general, systematic effects do not act all in the same di- 
rection and require a very detailed modelling in order to cor- 
rect for an amount that might be not very significant. Dif- 
ferent effects can combine to give a nearly null bias, as sug- 
gested from analysis of subsamples. In particular, cooling 
flow clusters in our sample are equally split between those 
elongated and compressed along the LOS, suggesting that in- 
cluding a temperature profile or other effects would not dras- 
tically change the results on their average structure. 

Departures from an isothermal ICM in hydrostatic equilib- 
rium should not introduce serious errors. The reliability of 
isothermal models to predict expected X-ray or SZE observa- 
tions was tested against high -resolution, hydrodynamic clus- 
ter simulations (Flores et al. 2005). As long as one focuses 
on cluster regions that are less sensitive to recent mergers, as- 
suming isothermal gas allows accurate X-ray estimations even 
if the gas has a strong temperature gradient. SZE decrement 
maps can also be accou nted for althoug h not as successfully 
as for the X-ray case ( IFlores et alJ2005t . 

7. DISCUSSION AND CONCLUSIONS 

We have studied the intrinsic structure of clusters of galax- 
ies by combining X-ray and SZE observations and modelling 
them as ellipsoids of revolution. We find that clusters of 
galaxies cannot be described as a population of either purely 
oblate or purely prolate spheroids with a random inclination 
of the symmetry axis respect to the LOS. A mixed popula- 
tion is instead required to model the data, prolate-like shapes 
being slightly preferred. 

We have then compared our results with what is expected 
from a population of randomly oriented galaxy clusters, hav- 
ing the observed projected ellipticities of our sample. In con- 
trast with the intrinsic population expected under this hypoth- 
esis, we find a significant excess of clusters nearly aligned 
along the line of sight. Furthermore, high inclinations angles 
are observed mainly for oblate clusters, while most elongated 
clusters turn out to be only prolate. If we believe the hypothe- 
sis of randomly oriented polar axes to be correct, the behavior 
of the clusters in our sample could be explained by a combi- 
nation of two effects: a first one caused by selection effects, 
which favours the detection of clusters more elongated along 
the line of sight, and a second one due to the fact that an oblate 
or prolate ellipsoidal model might not be correct for galaxy 
clusters, and a more general triaxial morphology should in- 
stead be used. Even though clusters in our sample were orig- 
inally selected according to their luminosity with a selection 
threshold well above the detection limit, a selection bias on 
the basis of X-ray surface brightness can still persists for ex- 
tremely elongated clusters. The second effect would instead 
lead to more dramatic conclusions. 

Higher accuracy measurements and a better understanding 
of systematics are needed to confirm our results. Although 
approximating galaxy clusters as isothermal systems in hy- 
drostatic equilibrium with ellipsoidal matter distributions can 
be quite simplistic for some clusters, it can give a first in- 
sight on their intrinsic structure. Some n umerical simulations 
also seem to support this view (iFlores et al. 2005) . A nat- 
ural development could be using more sophisticated models 



than the isothermal-/? profile to describe the ICM. Due to the 
infinite extent assumed in this model, the slope of the sur- 
face brightness distribution at large radii could be too shallow 
and might miss a progressive steepening with radius. Unfor- 
tunately, data of sufficient quality are often missing too and 
ICM can be traced with sufficient accuracy only up to few core 
radii. Furthermore, the /3-model provide a very good frame- 
work for the gas density distribution of most of the observed 
clusters in our sample and still provide a bench-mark for com- 
parison in a statistical sample. High-resolution spatially re- 
solved spectral measurements of temperature profile are re- 
stricted within the very inner regions for most of the clusters 
and some extrapolation technique, with a degree of unavoid- 
able liberty, is therefore required (Schmidt et al. 2004). Fur- 
thermore, from recent studies there appeared to be little dif- 
ference in the accuracy of the mass estimation from X-ray 
or SZE measurements whether using isothermal or general, 
numerica lly or observati onally motivated, temperature profile 
methods JHalknan et al.l l2b05). While waiting for more pre- 
cise measurements on a cluster-by-cluster basis, using only a 
few numbers, i.e. a single temperature, a constant axial ratio 
and a simple model for the ICM distribution, is still a con- 
servative approach. A crucial step to strengthen the statistical 
significance of our analysis would be the use of an unbiased 
and larger sample. The sample we have been considering in 
this paper is restricted to X-ray selected clusters for which a 
SZE analysis had been already reported. Preferential inclu- 
sion of high-S/N clusters could bias the sample towards clus- 
ters highly elongated towards the LOS. 

The problem of determining the intrinsic shape of galaxy 
clusters could be in principle completely solved for quite 
general morphologies, but uncertainties on observed quan- 
tities make this task still hard to obtain. An improved 
accuracy in the measurement of the SZE would also al- 
low to look for trends in the intrinsic structure with mass 
or redshift as shown by numerical simulations (Jing & Sut^ 
2002; Kasu n & Evrar d 2005; Paz et al. 2006). Unfortunately, 
present data are n ot sufficiently precise to test such predicted 
co rrelation s jDe Filippis et al. 2005J). As we have discussed 
in Paper 1", a proper triaxial structure is well suited to fit ob- 
servations. Additional independent constraints from gravita- 
tional lensing observations can break the degeneracy on the 
intrinsic shape of galaxy clusters and discriminate between 
triaxial spheroids and ellipsoids of revolution. In view of 
future progresses on the observational side, a more accurate 
adaptation of our method, in which both the assumptions of 
isothermality and of a density profile approximated by a f3- 
model are abandoned, will be presented in a forthcoming pa- 
per. Our complete 3-D deprojection method will also be in- 
vestigated with numerical simulations. The comparison of 
the projected properties of 3-D simulated haloes with obser- 
vational data could help to understand the intrinsic structure 
and to test the approximations and hypotheses used JPaz et 
l2006h . With respect to similar previous analyses of groups 
based on the optical distribution of member galaxies, our 
method would have the advant age to be not af fected by sys- 
tematics due to finite sampUng JPaz et al.l2006h . 
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APPENDIX 

TWO-DIMENSIONAL PROJECTION OF A THREE-DIMENSIONAL ELLIPSOID 

The expressions reported in Sections|2]and|5lcan be obtained from the results in lPaper j . Let us denote the observer coordinate 
system as {xi ohs}, i — 1,2,3. The polar axis forms an angle i with the line of sight, x^^ohs- We assume that the major axis of the 
projected ellipses always lies along the xi.obs-axis, i.e. the polar axis lies in the xi^obs — 2:^3, obs plane in the prolate case and in 
the X2,obs — 2:3 obs plane in the oblate case. As an example, in the prolate case, the relation between the two coordinate systems 
are 



a;i,int : 

2^2, int ^ 
a;3,int ^ 



= a;i,obs COS I 

= a^2,obs, 

= 2:3, obs COS i 



a:^3,obs smz 



obs smi 



Following the notation in lPaper 3. for a prolate ellipsoid vi = V2 = e-mt 
dc3 = ^c/eint. The observed core radius is 

( ^ 

= erx 



(Al) 
(A2) 
(A3) 

^'c3 = dc', for an oblate ellipsoid wi = U2 = l/eint, 



'proj 



Cint 



prolate case 



1 oblate case 

Substituting in the formulae reported in lPaper j we obtain the expressions used in Sections|2and|5] 



(A4) 
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